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Abstract 

We study quantum communication through an anisotropic Heisenberg XY chain in a 
transverse magnetic field. We find that for some time t and anisotropy parameter 7, 
one can transfer a state with a relatively high fidelity. In the strong-field regime, the 
anisotropy does not significantly affect the fidelity while in the weak-field regime the 
affect is quite pronounced. The most interesting case is the the intermediate regime 
where the oscillation of the fidelity with time is low and the high-fidelity peaks 
are relatively broad. This would, in principle, allow for quantum communication in 
realistic circumstances. Moreover, we calculate the purity, or tangle, as a measure 
of the entanglement between one spin and all the other spins in the chain and find 
that the stronger the anisotropy and exchange interaction, the more entanglement 
will be generated for a given time. 
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1 Introduction 



The transfer of a quantum state from one place to another is an important 
task in quantum information processing. A quantum state prepared by one 
party needs to be measured by another party at a distance. Long distance 
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communication between two parties, for example, in quantum key distribu- 
tion p], can be realized by means of photons. In this case, photons have the 
advantage that they have an extremely small interaction with the environ- 
ment and also travel long distances quickly through optical fibers or empty 
space. However, for short distance communication, such as connecting distinct 
quantum processors or registers inside a quantum computer [2113] . conditions 
and requirements are different [I] . Recently, quantum communication through 
spin chains has been intensively investigated for this purpose [5|6lT7] . 

The primary scheme is that one quantum state is produced at one end of the 
chain; it evolves naturally under spin chain dynamics; and at some time t we 
receive the state at the other end [5] . For more complex systems, Christandl et 
al. suggested a perfect state transfer algorithm which can transfer an arbitrary 
quantum state between two ends of a spin chain JB], or a more complex spin 
network [7]. In addition, researchers have investigated measurement-assisted 
optimal quantum communication by a single chain [8] and parallel chain [9|[T0] . 
the entanglement transfer through a Heisenberg XY chain [TT|T2] and parallel 
spin chains [13J, enhancement of state transfer with energy current [14] and 
entanglement transfer by phase control [ToTTB] . However, due to the complex- 
ity of the problem researchers usually consider cases where the magnetization 
(z component of the total spin J2i Sf) is a conserved quantity, which means 
that Sf, H] = 0. For the cases that the Hamiltonian does not satisfy above 
conditions, L. Amico et al. studied the dynamics of entanglement and found 
that the anisotropy of the Hamiltonian has an evident effect on the evolution 
of entanglement [TT1[T8] . They calculated the entanglement using the out-of- 
equilibrium correlation functions, but did not investigate state transfer in these 
systems. The anisotropy and magnetic field effects on the entanglement trans- 
fer in two parallel Heisenberg spin chains was later investigated in Ref. [19] . For 
the simple cases where each chain only has two spins, it was determined that 
perfect entanglement transfer can be realized between spin pairs by adjusting 
the magnetic field strength and the anisotropy parameter. 

In this paper, we study quantum communication in an anisotropic Heisenberg 
XY chain with small number of particles N and an arbitrary initial state. 
Following the scenario of earlier work on state transfer through spin chains, 
we encode a state to be transferred at the first site of the chain and, without 
external control, let the chain freely evolve. After time t the state is to be 
readout at the r th site of the chain. But different from the other cases with 
\J2i Sf, H] = 0, the quantum communication channel, the spin chain, not only 
transfers the state, but also generates entanglement. 

This paper is organized as follows: In part two we will give our model Hamil- 
tonian and calculate the time dependence of the one-site correlation function 
which will be used in the expression of the fidelity and tangle. In part three, 
we discuss the time evolution of transmission fidelity and purity in a short 
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chain for three regimes: strong-field, weak-field and intermediate regime. The 
final part is devoted to the conclusions. 



2 The model and the calculation 



The Hamiltonian of the anisotropic Heisenberg XY chain in a uniform trans- 
verse magnetic field h is given by [20] 



N N 

# = _ ./".s/-.s/. ; + j y s y s y +1 ) - J2 ns;. (i) 
i=i i=i 

where S" is the a spin operator (a = x,y,z) at site i, J x and J y are the 
anisotropic exchange interaction constants, and h is the transverse magnetic 
field. We assume periodic boundary conditions, so that the N th site is identified 
with the th site. The standard procedure used to solve Eq. ([1]) is to transform 
the spin operators S°" into fermionic operators via the Jordan- Wigner(J-W) 
transformation 

a = 5r,c„ = (-2S?)(-2£|) ■ ■ ■ (-23U)S-,n= 2, ...N, (2) 
where c n are one-dimensional spinless fermions annihilation operators and 

It will be convenient to introduce operators Ai = c\ + q, and B\ = c\ — q, 
which fulfill the anti-commutation relations 

{A l ,A m } = -{B l ,B m } = 2S lm , {A h B m } = 0, (3) 

In terms of these operators, the J-W transformation reads 

S? = \a!^\A s B 8 , Sf = - l -B l l l[A s B s , 5f = ~^S,, (4) 

A s=l A s=l A 

Using the J-W transformation, Eq. ([1]) becomes the bi-linear form 
H = - [( c Ui+i - Ci4+i) + 7 (44+i - QQ+i)] + h (c\ci - 2) } ■ ( 5 ) 

where J = |( J x + J y ), and 7 = (J x — J y )/(J X +J y ) is the anisotropy parameter. 
The limiting values, 7 = and 1 correspond to the isotropic and Ising chain, 
respectively. 

Eq. §5§ can be diagonalized by the transformation [TT]n%] 

Vk = iE e ' H (^ c i + Who]), (6) 
viV , 
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where 



1 — (h + J cos k)/X k 



/3k = sign(J , y sin k) 1 



'l + (h+ Jcos&;)/Afc 



2 

A fc = v /(/i + J cos A;) 2 + JV sin 2 /c, (7) 
The Hamiltonian in Eq. ([5]) then becomes 



H = 12 Xk (vivk - ^ 



where the wave number k = 2irm/N with —N/2 < m < N/2. 

Now that the XY Hamiltonian has been diagonalized, we will calculate the 
evolution of the operators c],Cj, which will be used to obtain the final state. 
From Eq. ([6]) and its inverse 

= 4^E e_i *W(*) - iforiUt)), (9) 
VN k 

where r)k(t) = e~ lXkt f]k(0). 
Then from Eqs. (JHD and © 



c j (t) = Y,[ai J (t)c l + b lj (t)cj} 



i 

J(+\ — i ~t m Ji 



cj(*)=EW^ + Wl. ( 10 ) 



where 



= ^ E e^'V A *< - 2za 2 sin A fc t], 
(*) = 4 12 e ik0 ~ j) <*kPk sin A fe t. (11) 

k 

For 7 = 0, and = 0, the evolution of the creation operator is 



/ 



In this case the ^-component of the total spin S z commutes with the Hamil- 
tonian and is a conserved quantity. 



4 



Now we seek to transfer a quantum state from the site s to r. First we assume 
that all the spins of the system are initially in spin down states. Then we 
encode the state |</?(0)) = a |0) + /3 |1) at the first spin of the chain. The initial 
state of the whole system is then |$(0)) = (a + (3c\) |0), where |0) denotes 
the state with all the spins down. For simplicity, a and (3 are taken to be real 
with a 2 + (3 2 = 1. Now our task is to calculate the fidelity of transmission of 
a state from the first spin to the r th spin of the chain. 

The reduced density matrix of the r th spin of the chain can be constructed 
using the operator expansion for the density matrix of a system on N spin- 1/2 
particles in terms of tensor products of Pauli matrices. The reduced density 
matrix is [21] 

/ \ + (s*) {s x r ) - % (sy) 



Pr(t) 

where (S?> y ' z ) means ($(0)| Sp v ' z |$(0)) . 



(12) 



The fidelity between the received state p r {t) and the initial state <p{0) is defined 



by F = J(<p(0)\p r (t) \<p(0)) which is 




-aP)(Sf)+2a/3(S»). (13) 

So if the parameters a = (3 = l/y/2, the fidelity F = \J\ + (Sf), which only 
depends on the value of (S®). 

Another quantity we want to calculate is the purity (also known as the tangle 
or one-tangle), which provides a measure of the entanglement between the 
spin at one site and the rest of sites in the chain. The purity is often expressed 
as 1 — Tr(p 2 ), but can also be expressed as 

r[p (1) ] =4det[p«], (14) 

where pW is the one-site reduced density matrix, Eq. (fl2l) . This is referred to 
as the one-tangle [18], which was apparently motivated by the tangle defined 
in Ref . [22] . Note that Eq. (fT4"I) is gives a valid measure of entanglement when 
the whole system is in a pure state and given there is only one parameter 
for a any such measure for a two-state system, this is as good as any other. 
However, if the system is in a mixed state, one must use some other measure 
of entanglement. The one-tangle, or purity, is connected to the Von Neumann 
entropy of the reduced density matrix through the relation 

% (1) ] = /4(1 + Jl-T\pW])), (15) 
where h(x) = —x log 2 x — (1 — x) log 2 (l — x). From Eq. ( TT2l) . the tangle/purity 
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can be written as 



r[p«] = 4det[p«] = 1 - 4«^> 2 + (S^) 2 + (S z r ) 2 ), (16) 

where the components S a are the components of the Bloch vector. 

Now we will calculate (S®(t)) (a = x,y,z) for use in Eqs. ([TBI and (IT5|) . 

(S?(f)> = « 2 (0| S?(t) |0) + /3 2 (1| S?(f) |1) + a/3«0| S x r (t) |1) + (1| S?(t) |0». 

(17) 

However, note that (0| S?(t) |0) = (1| S?(t) |1) = 0. 
Using Wick's theorem, 



(S?(t)) = i«/3{(0| (A lC \ + Cl A x ) |0) (0| B 1 A 2 S 2 ...A._ 1 S r _ 1 A. |0) - 

(0| (5x4 + Cl Bx) |0) (0| AiAa^.-A-i-Br-iA. |0> + 

• • • + (0| {A r c\ + c x Ar) |0> (0| A 1 B 1 A 2 B 2 ...A r _ 1 B r _ 1 10)}. (18) 



For the simple case, N = 3, r = 2 



<^(*)> = -f {(0| (A lC j + crAx) |0) (0| S X A 2 |0) 

-(oKs^ + dSOio) (o|aa 2 |o) 

+ (0| (A 2 cl + ciA 2 ) |0) (0| AiSi |0) . (19) 

The value of (Sf) can be obtained from the above expression by replacing 
A r -> £ r and 1/2 -»• -i/2. 

Now we need to calculate contractions of two field operators (0| AjB m |0), 
(0| A,A m |0), (0| Bj-B m |0), (0| B^An |0). From Eqs. (HDD- (ED, we get 
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(0\A j (t)B m (t)\0) = 5 jm 

4 

k 

+a k (3 k (l - 201) sin k(j - m)\ sin 2 \ k t 



2 

+ — ^2a k p k cos k(j - m)sin2X k t, (20) 

k 



{0\Mt)A m (t)\0) = 5„ 



jm 

Ai 

-If Et^Ai 1 ~ 2a l) cosk(j - m) 

k 

—2a\(3 k sin k(j — m)] sin 2 \ k t 
2i 
"N 



2i 

+ — ^a k f3 k sink(j -m)sin2X k t, (21) 

k 



and 



{Q\Bj{t)B m {t) \0) = -5 jm 

Ai 

-^EM^ 1 ~ 2a l) cosk(j - m) 

k 

+2a\0l sin k(j — m)] sin 2 X k t 
2i 

+ — ^2 a k (3 k sin k(j - m) sin 2\ k t. (22) 



(0| Aj(t)c\ + dA^t) |0> = - ^ [cos fc(l - j) cos A fc t 

— (1 — 2a 2 ) sin k(l — j) sin Afct 
+2a fc /5 fc cos/c(l - j) sinAfct], (23) 

(0| B 3 {t)c\ + Cl Bj(t) |0> = ^EK 1 - 2a 2 )cosA;(l - j) sinA fe t 

k 

+ sin fc(l — j) cos Afci 

+2a fc /5 fc sin fc(l - j) sin \ k t] . (24) 



Now we can calculate (S*) 

( s r) = -\\ a2 (0| A-^r |0) +/3 2 (1| A r B r \l)+ap((0\ A r B r |1) + (1| A r 5 r |0))]. 

(25) 

Note that, since (0| A r B r |1) = (1| A r B r |0) = 0, then (0| A r (t)B r (t) |1) = 



[l\A r (t)B r (t) |0) = pm 



(0\A r B r \0) + P A |6i r | -|ai r | 



|2 



(26) 



Now, from Eq. (120]) . (0| A r S r |0) = 1 - £ £* sin 2 A fe t. 
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For the case 7 = 0, = 0, we find that 



a(5 



£>s[A;(r-l)-A fc t], (S z r ) = p 2 a lr 



k 



2 



1 



(27) 



which agrees with Bose's case |5]. The a\ r gives the transmission amplitude 
fr 8 {t) of an excitation (the |1) state) from the 1 st to the r th spin. When 
a — j3 — 1/ a/2, the fidelity is 



which agrees with Ref. [14] . 



3 Results and discussion: Fidelity of state transfer 

We will now investigate the performance of a spin chain for which there exists 
some anisotropy parameter 7 and the ^-component of the total spin is not a 
conserved quantity. We first perform numerical calculations of the fidelity F 
for a particular, and small N value and discuss the variations with changes in 
physical properties of the system. 

As stated above, after some time t, the state at the r th spin will be measured, 
where r = N/2 + 1 for N even and r = (N + l)/2 + 1 for N odd), i.e. the 
sending and receiving positions are at opposite sites on the chain with periodic 
boundary conditions. Clearly the complexity of the computation grows with 
the value of r because the expansion of (S*(t)) has (2r — 1)!! terms. So for 
simplicity, we take N = 5, r = 3 as an example. We will also choose a particular 
state (y/3/2 |0) + 1/2 |1)) to analyze. However, it is important to note that we 
have examined several states and have found these trends typical; they exhibit 
similar variations, only the maximal values of the fidelity and purity/tangle are 
different. In each case the final state Eq. (fl~2l has a time dependence described 
by Eqs. ffTTj) and fl25|) . The state to be transferred is encoded at the first site of 
the chain. Then as time evolves the state propagates to the r th site. The time 
evolution of the operators (S°(t)) (a = x, y, z) can be expressed as a product of 
all contractions of the two operators A(t) and B(t), which are superpositions 
of many different states in the basis spanned by the Jordan- Wigner fermions 
with definite momenta [23] . This is clearly more complex than the isotropic 
case (7 = 0) where the state transfer can be characterized in terms of their 
dispersion [5]. 



(28) 
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Fig. 1. The fidelity F is plotted versus time t and anisotropic parameter 7 for (a) 
Strong-field regime, h = 1.0, J = 0.1. (b) Weak-field regime, h = 0.1, J = 1.0. (c) 
Intermediate case, h = J = 0.5. The number of sites N is 5, the input state is 
encoded in the first site and output state is at the 3 rd site, a = y3/2,/3 = 1/2. 



3.1 The strong-field regime 



The fidelity F as a function of time t and the anisotropic parameter 7 are 
shown in Fig. [T] with different parameter values for J and h. First in the 
strong-field regime (Fig. [H(a)), J/h <C 1, the dispersion is given by = 
h[l + (J/h) cos A; + (J//i) 2 (cos 2 k + 7 2 sin 2 k)/2 + o((J/hf)\. The effect of the 
anisotropy can be neglected if we only consider the term which is first order 

in J/h. The parameters in Eq. (J7]) become — > 0,/?/. — > ±1 and bij — > 
in Eq. ffTTl) . So the contractions of two field operators in Eqs. fl2D]) - (1221) are 
negligible for j 7^ m which indicates that using only a uniform magnetic field 
(global interaction) cannot cause a uncorrelated state to become correlated. 
In Fig. [T^a), it is clear that the anisotropy does not significantly affect the 
fidelity, but the presence of the cosine term in the first order of the dispersion 
produces the observed oscillation of fidelity with time t. 
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3.2 The weak- field regime 



In the weak-field regime /i/J<l (Fig. QJb)), the dispersion can be written as 
Xk = J[2(h/J) cos k + (h/ J) 2 + cos 2 k+"f 2 sm 2 k] 1 ^ 2 . In this case, the anisotropy 
has pronounced effects on the fidelity. We see that increasing the anisotropy 
does not always decrease the fidelity. For certain values of the parameters 7 
and t, the fidelity is greater, which corresponds to a constructive interference. 
For example, there are peaks with 7 7^ such as (7 = 0.28, t = 27.70, F = 
0.98) and (7 = 0.42, t = 7.10, F = 0.98) that exhibit this behavior. With 
increasing 7 and t, the frequency of oscillation of F becomes greater. The 
higher-frequency oscillation can be attributed to the fact that with increasing 
7 the initial state differs more from the true ground state, and therefore, with 
the evolution of time, it exhibits fluctuations, even near the ground state. 



3.3 The intermediate regime 



The intermediate regime (J ~ h), which is the most interesting, is shown in 
Fig. d](c). For relatively short times < t < 10, the anisotropy does not have a 
pronounced effect on the fidelity. The oscillation of F with time t is relatively 
slow in this case compared to the cases in both the strong- and weak-field 
regimes, Fig. [T](a) and (b) respectively. When the strength of the exchange 
interaction J or the magnetic field h dominate, the oscillation of F with time 
t is greater. But for intermediate regimes, the competition between J and 
h gives a smaller oscillation frequency for F. Furthermore, the high-fidelity 
peaks are fairly broad and therefore correspond to values of the fidelity which 
are stable under perturbation of the parameters. 



3.4 Comparison 



In summary, realistic quantum communication devices will require larger fi- 
delities in a shorter times. One may well have expected a significant oscillation 
of fidelity with time and anisotropy in the intermediate regime. However, this 
is not the case, but rather the intermediate regime has the highest fidelities 
in the shortest amount of time with the more stable values for the fidelity. So 
if one wants to achieve the highest fidelity under our assumed conditions, the 
intermediate regime is, surprisingly, the best choice. 
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Fig. 2. The one-tangle, or purity, at the 3 site as a function of time t and 7. 
The initial state of the system is (y/3/2 |0) + 1/2 |1)) <g> |0). (a) Strong-field regime, 
h = 1.0, J = 0.1. (b) Weak-field regime, h = 0.1, J = 1.0. (c) Intermediate regime, 
h = J = 0.5. 

4 Results and discussion: Entanglement 

It is clear that the initial state of the system is an unentangled state. However, 
with the evolution of time, entanglement is generated in the chain. In order to 
quantify this change, we have calculated the purity, or tangle, at the 3 rd site 
which measures the entanglement between the 3 rd site and all the other sites 
in the chain. The tangle at 3 rd site as a function of time t and anisotropy 7 
is plotted in Fig. [2] and Fig. [3] with different initial states (a |0) + (3 |1)) ® |0) 
and vacuum state, respectively. From Fig. [2], for small anisotropy (7 < 0.1), 
the tangle is negligibly small. In the strong- field regime, from Fig. [3]^a), the 
anisotropy does not have a significant affect on the tangle similar to behaviour 
of the fidelity as seen in Fig. [21(a)- Also, the strong-field tangle is relatively 
small compared to the weak-field regime and the oscillation of tangle with 
time t is suppressed. Furthermore, the time for the tangle reach its first peak 
in the strong-field regime is at t « 45 while for the intermediate regime it is 
t » 5. 

In the weak-field regime, which is plotted in Fig.[2]^b) and Fig.[3^b), we see that 
a stronger exchange interaction even with a stronger anisotropy can generate 
more entanglement at certain times. For example, at 7 = 1.0 and t = 1.5, 
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T = 1.0. 



In the intermediate regime, the behaviour of the tangle is similar to the weak- 
field regime. With increasing anisotropy, the tangle increases. However, even 
in the 7 = 1.0 case, it does not reach its maximal value r = 1.0. Also, as 
J/h increases, the tangle oscillates more rapidly, but different from the effect 
of anisotropy on fidelity in the intermediate case, (see Fig. [2(b)) the tangle 
increases with increasing 7. This is an effect of entanglement dynamically 
generated from the ground state and depends on the anisotropy. When the 
initial state is the ground state, the system will still generate entanglement, 
which exists only in the 7 7^ case. This is due to the double spin-flip operator 
terms 7(cJcJ +1 -Qc m ) in Eq. p. In this case, (0| S*(t) |0) = (0| S^t) |0) = 0, 
and r[p«] = 1 — 4 <0j S*(t) |0) 2 . Unlike Fig. [2J the one-tangle always equals 
zero when 7 = since the vacuum state is the ground state of the system. 
And the r th spin will always remain in the spin-down state and will never 
be entangled with the other spins. For the two initial states of the system, 
Fig. [2(b) (c) and Fig. [3(b) (c), we find that the time evolution of the one- 
tangle shows similar behavior with increasing 7, which means the one-tangle 
is not sensitive to the initial state of the system when a strong anisotropy is 
present in the case of weak-field and intermediate regimes. This shows that 
the anisotropy parameter aides in the generation of entanglement in the spin 
chain. 



Fig. 3. Same as Fig. [2] except that the initial state is the vacuum state, |0). 
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5 Conclusions 



In conclusion, we have investigated quantum state transfer through an anisotropic 
Heisenberg XY model in a transverse field and also entanglement generation 
in that same system. The interest in these problems stems from the possible 
use of spin chains as a communication channels. We expect that a realistic 
ferromagnetic material would have some anisotropy and have shown that this 
anisotropy can have significant effects on the fidelity of state transfer as well 
as the entanglement in the chain. 

Specifically, we have calculated the fidelity and the one-tangle, or purity, for 
three different cases: a weak external magnetic field, an intermediate regime, 
and strong external magnetic field. We found that in all three cases a relatively 
high fidelity can be obtained for certain times and vales of the anisotropy. How- 
ever, in the intermediate regime, the oscillation of the fidelity with time is fairly 
low, and the peaks fairly broad. Furthermore, a fairly high fidelity is achieved 
in a relatively shorter time for the intermediate regime. This would imply 
a more reliable output, in a shorter time, when some anisotropy is present, 
compared to when there is none. Thus the intermediate regime presents some 
interesting and somewhat surprising results for state transfer and is also the 
best choice for reliable transfer. 

We have also calculated the one-tangle, or the purity of a typical one-particle 
state in the chain. We began with a pure initial state and found that the 
stronger the anisotropy and exchange interaction, the more entanglement will 
be generated. This indicates that anisotropy also aides in the production of 
entanglement in the chain. 
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